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Abstract — In the multiuser spatial multiplexing multiple-input 
multiple-output (MIMO) system, the joint transmitter-receiver 
(Tx-Rx) design is investigated to minimize the weighted sum 
power under the post-processing signal-to-interference-and-noise 
ratio (post-SINR) constraints for all subchannels. Firstly, we show 
that the uplink-downlink duality is equivalent to the Lagrangian 
duality in the optimization problems. Then, an iterative algorithm 
for the joint Tx-Rx design is proposed according to the above 
result. Simulation results show that the algorithm can not only 
satisfy the post-SINR constraints, but also easily adjust the 
power distribution among the users by changing the weights 
accordingly. So that the transmitting power to the edge users in 
a cell can be decreased effectively to alleviate the adjacent cell 
interference without performance penalty. 

Index Terms — spatial multiplexing, MIMO, power allocation, 
Lagrangian duality. 

I. Introduction 

Spatial multiplexing for the multiple-input multiple-output 
(MIMO) systems, employing multiple ttansmit and receive 
antennas, has been recognized as an effective way to improve 
the spectral efficiency of the wireless link [1|. More recently, 
the multiuser schemes have been investigated for the spatial 
multiplexing MIMO systems. This paper focuses on the down- 
link multiuser schemes in which each user can not cooperate 
with the others thus suffers from the interference from them. 

Mainly, there are two kinds of multiuser schemes. One is the 
precoder or the transmit beamforming, such as the dirty-paper 
coding (DPC) [2| and the zero-forcing (ZF) 0, etc., which 
mitigates the multiuser interference only by processing at the 
transmitter. The other is the joint transmitter-receiver (Tx-Rx) 
design, such as the nullspace-directed SVD (Nu-SVD) |4| and 
the minimum total mean squared error (TMMSE) [5|, etc. 
In general, the former possesses lower complexity but more 
performance penalty. With the great development of signal 
processors, the latter gradually draws more attention. 

For the joint Tx-Rx design, the schemes proposed in 
|4||5] minimize mean squared error (MMSE), or maximize 
the capacity under the transmit power constraint. Whereas 
on some occasions, such as the multimedia communication, 
it is required to minimize the total transmit power while 
guarantee the quality of service (QoS). ||6)0 investigate the 
beamforming and the power allocation policy when all users 
are subjected to a set of post-processing signal-to-interference- 
and-noise ratio (post-SINR) constrains in the uplink SIMO 
and the downlink MISO. [8|[9] extend this work to the 



downlink MIMO and the MIMO network, however the MIMO 
systems discussed in [8] [9] are assumed that there is only 
one substream between each pair of the transmitter and re- 
ceiver. In other words, only the multiuser interference appears 
in the so-called diversity MIMO system in [8|[9|. For the 
multiuser spatial multiplexing MIMO system, however, both 
the multiuser interference between individual users and self- 
interference between individual substreams of a user should 
be mitigated. 

For the downlink, the transmit beamforming affects the 
interference signature of all receivers, whereas the receive 
beamforming only affects that of the corresponding user. 1 7 1 1 8 ] 
construct a dual system, called the virtual uplink, and indicate 
that the virtual uplink can obtain the same post-SINR as the 
primary downlink. Moreover, the receive beamforming matrix 
of the virtual uplink is identical with the transmit beamforming 
manix of the primary downlink. The design of the downlink, 
therefore, can resort to the virtual uplink. 

In this paper, we extend the duality derived for MIMO 
network in J9) to the multiuser spatial multiplexing MIMO 
system. According to the uplink-downlink duality, we propose 
a joint Tx-Rx scheme to minimize the weighted sum power 
under the post-SINR consttaints of all the subchannels. 

Notation: Boldface upper-case letters denote matrices, and 
boldface lower-case letters denote column vectors, ir(-), (•)*, 
(■) H , || • lb an d || • | |f denote trace, conjugate, conjugate trans- 
position, Euclidian norm and Frobenius norm, respectively. 
diag(x) denotes a diagonal matrix with diagonal elements 
drawn from the vector x. denote the (i,j)-th element 

and j-th column of a matrix, respectively. 

II. System model 

We consider a base station (BS) with M antennas and 
K mobile stations (MS's) each having Ni(i = 1,...,K) 
antennas. There are Li(i = 1,...,K) substreams between 
BS and MS,(i = 1, . . . , K), that is to say, BS transmits L t 
symbols to MS; simultaneously. The signal recovered by MSfe 
can be written as 

A' 



y° k L 



Af H fc ^2 ^idiag(^/pl)xi 
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(1) 



where y® L 6 C LkXl is the recovered signal vector, e 
C LiXl (i = 1, . . . , K) is the transmitted signal vector from BS 
to MS; with zero-mean and normalized covariance matrix I. 



Pi € lZ LiXl denotes the power vector allocated to MS^. A lin- 
ear post-filter Afe £ C NkXLk is used to recover an estimation 
of the transmitted signal vector x/j . The MIMO channel from 
BS to MSfc is denoted as H fc e c NkXM , and assumed flat 
faded. Hence, its elements are the complex channel gains, and 
they are independently identically distributed (i.i.d.) zero-mean 
complex Gaussian random variables with the unity variance. 
Moreover, the perfect channel state information are assumed 
available at both transmitter and receiver via some way, for 
example, channel measurement at receiver and fast feedback 
to the transmitter for the frequency division duplex (FDD) 
systems, or invoking the channel reciprocity in time division 
duplex (TDD) systems. S C MxLi is used to weight x^ and 
transform it into a M x 1 vector, S C NkXl is the noise 
vector with the correlation matrix R„ = u^I. For simplicity, 
in the sequel we assume L\ = . . . = Lk = L. 

We design the A^, and Pk(k — 1,...,K) in (1) to 
minimize the weighted sum power under the post-SINR con- 
straints, which can be denoted as the following optimization 
problem. 



(4) 



j-th substream can be denote as 

„H ns,M„ 

qtatrDL _ th>2—M '° 
OtHUkj - I+n , DL 

a /.., I1 7.., a <*•••.; 
*IT DL = E PMH fc b M b^Hf+ 

K 

J2 H fc B m ^a 5 (p m )B^Hf + u 2 n \ 

If x m = [x m ,i,...,x m) i], yf L = [yk,i,---,yk,L], the link 
power gain between x m>n and y k j can be denoted as 



[ ( t ) k,j]m,n — | | a fe^jHfcb mjra | || 



then (4) can be rewritten into 

PkJ [<frkj]kj 



SINR°f 
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(6) 
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P,A fc ,B fc 

s.t. SINRgj > lK 



(k = l,...,K,j = l,...,L) 



(2) 



By substituting (6) into the constraint inequality of (2), we 
obtain 

cLp + ^Ha^lll^O (k = l,...,K,j = l,...,L) (7) 



where p = [pf , . . . , p L K ] 1 and w e 1Z 



iKLxl 



is the weight where the m-th element of c k , e 7e XLxl is 



vector, w affects the power distribution among users, and its 
value is determined by various factors, such as the positions 
of users in a cell and the interference environment of the 
neighboring cells. 7^ is the given post-SINR goal for the 
MSfc's j-th substream. 



III. The proof of uplink-downlink duality 



If A fc = [a feil ,...,a fei i], B fc = [b k>1 , . . ., b k ,L], Pk 
[p k ,i, ■ ■ -,Pk,L] , (1) can be rewritten into 



7fc,j 



m = (k — 1)L + j 

[0fej]r ? ], m _ ( r™]_i )I , m^{k-l)L + j 

(8) 

where rounds ^ to the nearest integer greater than or 
equal to 22 . Write (7) into the matrix form, we obtain 

Cp + d < (9) 

where C <E K KLxKL and d e K KLxl are 

C = [c M , . . ., Ci. L , . . ., ck.i, . . ., c K , L ] T 

d = ^[||a 1 , 1 |||,.. ) ||a ll z||l,..,||a J c, 1 ||l,..,||a J c, i |lS] T 

(10) 

So, (2) is equivalent to the following optimization problem 
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(11) 
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B i diag(v / p7)xi + Af n fc 



(3) 

The diagonal elements of the first part in the right-hand side 
(RHS) of (3) denote the useful signals, and the non-diagonal 
elements denote the self-interference. The medial and the last 
parts in the RHS of (3) denote the multiuser interference and 
the noise, respectively. Moreover, the post-SINR of the MS/c's 



Subsequently, to obtain the Lagrangian duality of (11) J5], we 
divide the solving process of (11) into two steps similar with 
[ 1 1 . First, assuming Ak and Bfc(fc = 1, . . . , K) are fixed, the 
Lagrangian function of (11) is 

L(p,X,n) = w T p + A T (Cp + d) fi T p (12) 

where A > 0, fi > are the Lagrangian multipliers associated 
with the inequality constraints. Then the Lagrangian duality of 
(11) is 

maxmin L(p,A, u) 

p (13) 
s.t. A > , £t > 



According to the Slater's condition, (11) is equivalent to (13). 

Since the gradient of the Lagrangian function (12) with respect 

to p vanishs at optimal points, we obtain w T — fi T — — A T C. 

Substituting it into (12), we obtain min L(p, A, /x) = d T A. 

p 

Moreover, as A > and fi > 0, (13) can be rewritten to 



max 
A 

s.t. 



d T A 

C T A4 
A > 



w > 



(14) 



Similar with (6)-(9), substitute (10) into (14), we obtain 

(k = l,...,K,j = l,...,L) 



max d T A 
s.t. SINRUj < j k ,j 



(15) 



where 



SINR™ = 



^r L = e 

K 



(16) 
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where A = [Af , . . . , A^]. Furthermore, SINR^j is the post- 
SINR of MSfc's j-th substream in the virtual uplink 



K 



y u k L 



(17) 
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(14) maximizes the weighted sum power under the maximum 
post-SINR constraints, however, it has no physical meaning 
||9l . But it can be shown that (14) is equivalent to the following 
optimization problem 



min d T A 

A rrr (k = l,...,K,j = 1, 

s.t. SINRVL> lk j v 



(18) 

Theorem 1: At the optimal point, the post-SINR constraints 
in (15) and (18) are active. And the solutions of (15) and (18) 
are identical. 

Proof: Without any loss of the generality, we assume 
SINRUL ^ Jk,j- From (16), we can find Xk,j contribute 
to the numerator of SINR^j and the denominator of 
SINR^Jm £ k,n ^ j). In' other words, SINR.%% is a 
monotone increasing function of Xk,j, while SINR^ n (m ^ 
k,n j) is a monotone decreasing function of Xk,j- So 
increasing Xk,j until SINR^ — "fk,j> we obtain a larger 
d T A without breaking any post-SINR constraint. Likewise, if 
SINB%*[ > 7fcj, decreasing X k ,j until SINR%% = 7 fej , a 
smaller d T A is obtained. As a result, the constraints of (15) 
and (18) become a linear equations C T A + w = 0, and its 
solution is A* = — (C T ) _1 w. 



Similar with Theorem 1, at the optimal point of (11) p* = 
CM. 

Summarize the above statement, we obtain the following 
conclusion. 

Theorem 2: In the downlink multiuser spatial multiplexing 
MIMO system, if the transmit and receive beamforming matri- 
ces are and A^(fc = 1, . . . , K), respectively, as long as the 
following conditions are satisfied, the downlink optimization 
problem (2) is equivalent to the virtual uplink optimization 
problem (18). 

1 ) In the virtual uplink, the transmit and receive beamform- 
ing matrices are A^ and B^(fc = 1, . . . , K), respectively. 

2) In the virtual uplink problem (18), the weight vector w 
is the noise power vector. 

3) In the virtual uplink problem (18), the noise power vector 
d is the weight vector. 

When Afc and Bfc(fc = 1, . . . , K) are not fixed, (18) is a 
joint optimization problem denoted as 

min d T A 

A,A fc ,B fc (k = l K j = 1 L) 

s.t. siNR^>jk,r ' ' ,J ' ' j 

(19) 

Theorem 3: If the noise power vector in the virtual uplink 
is the weight vector w, the joint optimization problem (2) 
is equivalent to (19). At the optimal point, the beamforming 
matrices of the virtual uplink and the primal downlink are 
common. 

Proof: Let (B*,A*,p*) (k = l,...,K) be the global 
minimum of (2). According to Theorem 2, (19) has the 
solution (A£, B£, A*) k = 1, . . . , K. Moreover, this solution 
is definitely the global minimum. Otherwise, a better solution 
of (2) would be found by applying Theorem 2 again. So 
the virtual uplink and the primal downlink have the common 
beamforming matrices. 

The weight vector w decides whether A^. and B^ are used 
to strengthen the useful signals or alleviate the interference to 
other users. When a user's weight turns higher, its transmit 
power will decrease. In this occasion, it benefits to apply the 
beamforming to increase the signal gain, as the interference to 
other users is much less important. On the other hand, once 
the weight gets lower, the beamformer should try to suppress 
interference to others (9)- 

To mitigate the adjacent cell interference, we can increase 
the weights of edge users in a cell, which would induce the 
declining of the transmit power from the BS to them. In order 
to hold the post-SINR under this circumstance, obviously, the 
beamforming matrices would be used to boost up the signal 
gain. 

IV. The joint Tx-Rx beamforming scheme 

It is rather difficult to solve the joint optimization problem 
(2) directly. However, it is easy to obtain Ak(k = 1, . . . , K) 
in the primal downlink, and so does Bfc(fc = 1, . . . , K) in the 
virtual uplink. Moreover, it is proved in the previous section 
that the primal downlink is equivalent to the virtual uplink, 
and they have the common beamforming matrices Afe and 
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Fig. 1. Block diagram of the joint Tx-Rx beamforming scheme 



Bfc(fc = 1, . . . , K). Therefore, we divide the solving process 
into four steps shown as Fig. Q] 

When p and Bfc(fc = l,...,K) are fixed, optimize 
A k (k = 1, . . . , K) to maximize SINR^f(k = 1, . . . , K, j = 
1,...,L). Then observing (4), it is a generalized Rayleigh 
quotient problem, and its solution is 



a fej - = a ktj /\\a ktj \\2 



(20) 



where £ max (X,Y) is the dominant generalized eigenvector 
of the matrix pair (X,Y). When A and Afe(fc = 1, . . . , K) are 
fixed, in the same way, Bfc(fc = 1, . . . , K) can be obtained by 



J fej7 



The proposed algorithm is summarized in the following 



5 Ml2 

,s,UL -r,I+n,UL\ 
, r < ; 



(21) 



Initialize B k (k — 1, . . . , K) and p l ' randomly. 

Set the noise vector of the virtual uplink 
to w . 

n = 

1 ) Update in the primal downlink. 

a) Calculate A<: n+1) (k = 1, . . . , K) from 

B ( k n \k = 1, . . . , K) and p (n) using (4) (20). 



b) Calculate C (n) from B k n) (k = 1 



and A 



(n + l) 



(k = l,...,K) using (5) (8) (10). 



c) Solve A (n) = -((C (n) ) T )- 1 w 
2 ) Update in the virtual uplink. 

a) Calculate B k n+1) (k = 1, . . . , K) from 



fn+l) 



(A; 



and AJf^ffc = 1 



■ ( " +1 V 

c) Solve p (n+1) = -(C (n+1) )- 1 d 
n = n + 1 



K) and A (n) using (16) (21) . 
■on. B[ n+1 \k = 1,...,K) 
K) using (5) (8) (10) . 



b) Calculate C (n+1) f rom B<™ +1) (k = 1, . . . , K) 



3 ) Repeat 1) and 2) until 

EIIA^-A^I 

k=l fc=l 

In the simulation, we set s = 0.0001 



E IIBl' 



(n) 



B 



(»+i)i 



f < e 



By iteration, (Aj c " +1) ,B^ +1) ,p(™+ 1 ))(A: = 1,...,K) con- 
verges to the optimal solution to the optimization problem (2). 
Once any element in p(' i+1 ) is negative, which indicates 
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Fig. 2. Total transmit power versus SINR goal 7, when K = 2,3,4 
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Fig. 3. Total transmit power versus SINR goal 7, when K = 5,6,7,8. 



not be attained, 7/.J should be decreased to relax the post- 
SINR constraints. When w = 1, the proposed algorithm is 
similar with the one in ifTTIl . 

V. Simulation result 

In this section, we assume that a BS with 8 antennas 
(M = 8) is communicating with K MS's each with 2 
antennas, [N\ = . . . = Nk = 2). Also we assume that 
the number of substreams of each MS is the same, equal to 
2,(Lx = . . . = Ljc = 2). QPSK is employed in the simulation 
and no forward error coding is considered. The post-SINR 
goals for all substreams are j(jk,j — l)- Additionally we 
assume MSi is an edge user in a cell, according to the previous 
section, a higher weight should be assigned to it to mitigate 
the adjacent cell interference. Thus, the weight vector is set to 
w = [w, w, 1, . . . , 1] T , where w is the weight corresponding 
to the two substreams of MSi and w > 1. 

Fig. 12131 plot the curves of the total transmit power EPfe j 

k,j _ 

versus the post-SINR goal 7, when w = 5. In Fig. |2] 
K = 2,3,4, the system configuration satisfies M > KL, 
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Fig. 4. When K = 4, 7 = lOdB, the transmit power of MSi and total 
power versus the weight w. 

the multiuser interference, thus, can be effectively suppressed 
through the beamforming [4|. Under the circumstance, increas- 
ing the transmit power of any user has nearly no effect to 
the post-SINR of other users. Therefore, all the substreams 
can attain relatively high post-SINR. In Fig. [3] K = 5, 6, 7, 8, 
M > KL does not hold any more. Consequently, the multiuser 
interference can not be effectively mitigated, which means any 
enhancement in the transmit power of any user is very likely 
to deteriorate the post-SINR of other users. As shown in Fig. 
[3] with the user number increasing, the available post-SINR 
of each user is decreased. When K = 8, only dB post- 
SINR can be attained. In these two figures, the total transmit 
power increases with the number of users and the post-SINR 
goal 7. Especially when K = 7,8 and 7 > dB in Fig. 
[3] due to the residual multiuser interference, the slopes of the 
curves are much steeper than that in Fig. [2] where the multiuser 
interference is negligible. And the steeper the curves are, the 
more power would be paid for the unit increase of the post- 
SINR of each user. 

Fig. [4] shows the curves of the transmit power Y^Pij °f 

3 

MS 1 and the total transmit power J^Pkj versus the weight 

k,j 

w, when K = 4 and 7 = lOdB. The left vertical axis is 
corresponding to the transmit power of MSi and the right one 
is to the total transmit power. Obviously, as the w is increasing, 
the transmit power of MSi is decreasing while the total power 
is increasing, because the optimization object is to minimize 
the weighted sum power J2 w k.jPk,j- Moreover, when w 

k,j 

changing from 1 to 20, the transmit power of MSi decreases 
almost lOdB, however the total power increases only about 
1 dB, which demonstrates that the proposed algorithm adapts 
the power allocation policy very effectively with negligible 
penalty on performance. 

VI. CONCLUSION 

In this paper, we investigate the joint Tx-Rx design for 
the downlink multiuser spatial multiplexing MIMO system. 
We show, first, the uplink-downlink duality has the following 



characteristics: 1) In both of the primal downlink and the 
virtual uplink, the substreams can attain the same post-SINR 
goal; 2) The beamforming matrices are common in both of 
the primal downlink and the virtual uplink. 

Based on the duality, a joint Tx-Rx beamforming scheme is 
proposed. Simulation results demonstrate that the scheme can 
not only satisfy the post-SINR constraints which guarantee 
the performance of the communication links, but also easily 
adjust the power distribution among users by changing the 
weights correspondingly, which can be used to diminish the 
power of the edge users in a cell to alleviate the adjacent cell 
interference. 
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